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ABSTRACT 



We analytically investigate the modes of gravity- induced star forma- 
tion possible in idealized finite molecular clouds where global collapse 
competes against both local Jeans instabilities and discontinuous edge 
instabilities. We examine these timescales for collapse in spheres, discs, 
and cylinders, with emphasis on the structure, size, and degree of in- 
ternal perturbations required in order for local collapse to occur before 
global collapse. We find that internal, local collapse is more effective 
for the lower dimensional objects. Spheres and discs, if unsupported 
against global collapse, must either contain strong perturbations or must 
be unrealistically large in order for small density perturbations to col- 
lapse significantly faster than the entire cloud. We find, on the other 
hand, that filamentary geometry is the most favorable situation for the 
smallest perturbations to grow before global collapse overwhelms them 
and that filaments containing only a few Jeans masses and weak den- 
sity perturbations can readily fragment. These idealized solutions are 
compared with simulations of star-forming regions in an attempt to de- 
lineate the role of global, local, and edge instabilities in determining the 
fragmentation properties of molecular clouds. The combined results are 
also discussed in the context of recent observations of Galactic molecular 
clouds. 

Subject headings: ISM: clouds, ISM: structure, stars: formation 
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1. INTRODUCTION 



Stars primarily form in clustered environments inside of mol ecular clouds 



conta ining tens of thousands of solar masses of material or more ( iLada fc Lada 
20031 ). These star- forming molecular clouds must fragment to give rise to the 
observed initial mass function (IMF) of st ars, with a characteristic mass on the order 



of a solar mass (e.g. 



Bonnell et al 



20071 ). and cannot be solely dominated by large 
scale, global collapse modes. Gravitational collapse within a molecular cloud leads 
to the formation of individual stars; however, the mechanism driving the formation 
and shape of the IMF is poorly constrained. While turbulence is believed to play a 



key role in producing significant i nhomogeneity 



gravity must also contribute (e.g. 



Bonnell et al 



rom w hich fragmentation can occur, 



20071 ). In fact, the idea that gravity 



is driving supers onic motions (or "turbulence" ) in molecular c 



for a while (e.g., 



Goldreich fc Kwan 



1974 : 



Ferrini et a 



19861 ) . Beyond the monolithic collapse envisi oned by 



1983 



ouds has been around 



Falgarone & 



in their critique of 



Goldreich fc Kwan 



Zuckerman fc Evand ( 119741 ) 



Puget 



( 119741 ). recent numerical and analytical 



studies show that global gravity can indeed dri ve "turbulent" motion s without 



leading directly to "catastrophic" co. 



2002 



2008 



Ballesteros-Paredes et al. 



Heitsch fc Hartmann 



2007 



l apse (e.g., 



Hartmann et al. 



Vazquez-Semadeni et al. 



2001 



2007 



Hartmann 



Field et al. 



20081 ). Depending on the morphology of the cloud, the 
fragmentation that occurs during such a large scale collapse might also lead to an 
IMF distribution of stars. In this scenario, the observed supersonic motions of the 
gas are determined by the geometry of the cloud, which depends strongly upon the 
initial conditions of the cloud. 
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Spatially resolved observations of star-forming molecular clouds usually 
exhibit rather complex geometries, including substructu re often associated with 



lower dimensionality, such as sheets anc 



1979 



Bally et al. 



1987 



*ada et al. 



1999 



ilaments (e.g., 



Hartmann 



Schneider &: Elmegreen 



2002 



Johnstone et al. 



2003 



Lada et al 



2007 



Myers 



2009 



Molinari et alJboiot kndre et al.lboioh . Filamentary 



structures have also been noted to form in molecular cloud simulations, although 



their formation mechanism varies wi th the type o: 



decay ing supersonic turbulence (e.g. 



Padoan et al. 



model used, be it driven or 



2007 



Bate 



Vazquez-Semadeni et al 



Burkert fe Hartmann 



2004 



20091 ). flow collisions 



201 01), or global gravita tional accelerations (e.g., 



Hartmann fe Burkert 



20071 ) 



In the absence o 
spheres (e.g. 
sheet s (e.g., 



(e.g.. 



Larson 



Larson 



Stodolkiewicz 



Nakamura et al 



tur 



j ulent motions , the collap se properties of isothermal 



1969 



1985 



Penston 



1963 



1993 



Myers 



1969 



Ostrike 



Shu 



19771 ). infinitely large, isothermal 



2009 ), and infinitely long, isothermal cyli nders 



1964; 



Fiege &: Pudrit 



Inutsuka fc Miyama 



1992 



1997 



20001 ) have been thoroughly studied. 



Equilibrium configurations of infinite sheets perturbed to form successive, parallel 
filaments and infinite, isoth ermal filaments perturbed to form strings of d ense 



cores have also been found (ISchmid-Burgk 



1967 



Curry 



2000 



Myers 



20091). Yet, 



the collapse of finitely long cylinders and discs, and the subsequent interplay 
between global and local collapse modes, are much less well understood. Previous 
simulations of finite cylinders show that these cylinde rs can fragment into multiple 



co ndensations (e.g., 



of 



Bastion 



Burkert Sz Hartmann 



1983 



Bastien et al. 



19911 ). while the disc simulations 



(120041 ) show no sign of local collapse modes. These disc 



simulations are dominated by the global collapse of the disc and the subsequent 



-5 - 



build up of mass along the outer ring of the disc (see also lHsu et al.ll2010l ). 



In this paper, we compare analytic estimates for the global collapse timescales 
to the local collapse timescales of small density perturbations embedded in different 
dimensional objects. We analyze whether small density perturbations are capable 
of collapsing significantly faster than the entire cloud for spherical, disc shaped, 
and filamentary clouds to determine whether strong perturbations, possibly induced 
by turbulent motions, are required to fragment these differently shaped molecular 
clouds. We present our analytic results in §0 In £j3j we discuss these results in 
context of recent observations and simulations and we present a summary of our 
findings in §HJ 



MODES OF COLLAPSE - ANALYTIC THEORY 



To understand the importance of global gravity in molecular clouds, we consider 
the large scale gravitational forces at play in comparison to the growth of local 
perturbations within a cloud. For a u niform spher e of arbitrary radius, the collapse 



is quickest for the largest scale (e.g., 



Larson 



19851 ). since gravity is a long-range 



force and global collapse is inevitable provided the region contains more than a 
Jeans mass of material and support from effects such as magnetic fields and rotation 
are negligible. In two-dimensions, however, an infinite sheet can be constructed 
for which there is no global gravitational force in the plane and the collapse is 



dominated by objects with radii similar to the thickness of the sheet (IL arson 



A similar argument can be made for infinite one-di mensional cylinders ( 



1985 



Inutsuka fc Miyama 



1992 



19971 ). As noted by 



1985|) 



Larson 



Burkert &: Hartmann 



(120041 ) 
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the lack of global gravity in these systems depends entirely on the infinite nature of 
the sheet or cylinder. In the following we show how the global gravitational collapse 
depends on both the dimensionality and the finite nature of simple geometries: 
spheres, discs, and cylinders. 



2.1. Thermal Support 



Thermal support is always present in a molecular cloud and will oppose 
gravitational collapse. This thermal support is most effective on small scales and 
can be overcome on larger scales. The Jeans length is the size scale above which, in 
three dimensions, thermal pressure is no longer capable of supporting an isothermal 
cloud. The Jeans length, Rj, is given by: 



R, 



7T 

~G~p 



where c s is the isothermal sound speed, G is the gravitational constant, and p is 
the mass density of the gas. For isothermal, infinite sheets, a stable equilibrium 
configuration can be found for all surface densities such that thermal pressure 
counteracts gravity. In these equilibrium configurations, the scale he ight of the ga s 



is on the order of the Jeans length of the material in the midplane ( jSpitzer 



1942 



19781 ) and small perturbations in the plane of the disc w ill only collapse if they 



19851 ). For infinite, isothermal 



are larger than the scale height of the gas (iLarsonl 
cylinders, the characteristic radius of an equilibrium configuration and the size 
scale on which thermal support is effective are on the order of the Jeans length at 



the central axis of the cylinder ( IStodolkiewicz 



1963 



Qstriker 



196J). Unlike infinite 
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isothermal sheets, infinite isothermal cylinders have a critical mass p er unit length 



value above which no equilibrium configuration exists ( lOstriker 



19641 ). 



In the following sections, we only consider the effect of gravity and do not 
explicitly include any thermal effects in calculating collapse timescales, although we 
do interpret our results in light of the Jeans length being the minimum scale which 
can collapse. We also do not consider the effects of turbulence, rotation, magnetic 
fields, or deviations from isothermality. 



2.2. 3D - Sphere 

For a sphere of density p and radius i? 3> Rj, the instantaneous global collapse 
time, t g , is given by 

The above instantaneous global collapse time differs by less than a factor of two 
from the collapse timescale obtained by properly followi ng the changes in the 



20071 ). A 



strength of the gravitational acceleration over time (e.g., ICarroll fc Ostlid 
small region, Ri <C R, at the centre of the sphere with an enhanced density (1 + e)p, 
where e«l, might act as a seed for local collapse. Since the density is only slightly 
enhanced, the Jeans length for this region is very similar to the unperturbed Jeans 
length and the region is stable against local collapse when Ri < Rj. However, 
provided Ri > Rj, the local collapse time, ti, for this inner region approaches 

1/2 

(3) 



27rG(l + e)p 
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where we have ignored the residual support due to thermal pressure. Comparing 
these two timescales reveals the difficulty in fragmenting an almost uniform spherical 
molecular cloud: 



t, 



:i+e) 



-1/2 



t, 



(4) 



For small amplitude density perturbations, there is little time for the local collapse 
to proceed before the global collapse consumes the entire spherical cloud. In 
three dimensions, local collapse beats global collapse on ly if there are significant 



inhomogeneities (strong perturbations) in the medium. iTohlind (119801 ) showed 
that the growth of this perturbation is even slower, relative to the global collapse 
timescale, if the perturbation is induced after global collapse has begun. 



2.3. 2D - Disc 



2.3.1. Accelerations 



For an infinitely thin disc with constant surface density £ and radius R 3> Rj, 
the central acceleration at a distance of r > R from the center of this disc is 



''ext 



4GS 



k\*)-e(* 

r I \ r 



(5) 



where K and E are th e first and second complete elliptic integrals respectively 



flWvse fc Mavall 



194 2f) . The cen tral acceleration of a point inside of the disc, i.e. 



r < R, is ( IWvse fc Mavall 



1942h : 



l int 



AG E 



R 



(6) 
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For a perturbed disc with surface density (1 + e) E inside of Ri <C R, the net 
acceleration is 



°net — < 



4GE^ 



4GE- 



Mi) -Hi) 



+4 



+ e 



R 



K(-k)-Hw) 



r < Ri 



R x <r < R. 



(7) 



For the case of R\ < r < R, 



Burkert &: Hartmannl (12004 ) showed that these elliptic 



integrals can be expanded to lowest order to give: 



•-net 



vrGEr / eRR\ 



(8) 



Figure [T] compares the accelerations of a perturbed disc calculated from the 
above exact and approximate solutions. For simplicity, G, E, and Ri have all been 
set to 1, e has been set to 0.1, and R has been set to 100. The approximate solution 
of Equation M reproduces the exact accelerations of Equation [7] for intermediate 
radii between Ri and R, but underestimates the accelerations closer to Ri and R. 



2.3.2. Timescales 

The square of the instantaneous collapse time at r is 

2r 



2R 
ttGE 



1 + 



eRR\ 



(9) 



The global collapse timescale, t g , can be determined by considering the collapse 
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time when r = R. That is 



2R 



2\ - 1 



^ = + ■ < 10 > 

Given that the initial surface density perturbation is small, e C 1, and Ri <C R : 

eR\ 



R 2 



« 1, 



o 2x1 



Thus, the global timescale is not changed significantly due to the small perturbation. 
Equation [9] can now be re- written as 



2\ -i 



t 2 (r)= ft il + e -^) . (12) 



9 \ 1 r 3 

Note that the collapse time t(r) is equal to the global collapse time and independent 
of r unless e > 0. When there is a density enhancement, however, the local time for 
collapse is dependent on r and thus, smaller regions can collapse faster than the 
global collapse. 

Assuming that the perturbed region is large enough to be Jeans unstable, 
Ri > Rj, the timescale for collapse of a local region with radius R\ is approximately 



tf - n[l + -J . (13) 

For this to be significantly smaller than the global collapse timescale 

eR 

e » (14) 



- 11 - 



Thus, it is possible to find a combination of perturbation length and amplitude that 
results in significantly enhanced local collapse. For a small density perturbation, 
e = 0.1, the perturbation length must be approximately 100 times smaller than the 
disc's total radius in order for the local collapse timescale to be 3 times smaller than 
the global collapse timescale. Recognizing that thermal support in the disc requires 
-Ri > Rj for local collapse to occur, such a disc would need to have assembled more 
than 10 4 Jeans masses of material! Even for a large perturbation, e = 1, the disc 
would still need more than 100 Jeans masses for the local collapse mode to be 3 
times faster than the global collapse mode. 

Figure [2] shows the collapse timescales for a disc with Ri — G — £ = 1, 
e = 0.1, and R/Ri = 100. The solid line shows the timescales calculated from 
the exact accelerations while the dashed line shows the timescales calculated from 
the approximate accelerations given by Equation [HJ The timescales for collapse 
interior to a radius of R\ are approximately a factor of three faster than the collapse 
timescales for the outer parts of the disc and the timescale for collapse at R\, as 
calculated from the approximate accelerations, is roughly equal to the collapse 
timescales for most points inside of Ri as calculated from the exact accelerations. 

The exact accelerations lead to collapse timescales approaching zero at radii 
of R and Ri, while the approximate accelerations produce timescales that do not 
show this asymptotic behavior. These extremely small collapse times are artifacts 
of having an infinitely thin disc with infinitely sharp edges and no thermal support. 
Real discs have finite heights and smoother edges, which should suppress this 
asymptotic behavior. We will come back to the issue of these infinitely small 
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collapse timescales at the edges in Section [2[ 



Thermal pressure provides for support against collapse on scales smaller than 
the Jeans length, Rj, and thus, the above analysis suggests that the preferred size 
of the perturbed region is R\ ~ Rj, since this provides the greatest leverage on 
the local to global collapse timescale. This length scale is also associated with the 
scale height of a thermally supp orted disc and is close to the preferred scale for disc 



fragmentation in infinite sheets ( 1L arson 



19851 ). 



2.4. ID - Cylinder 

2.4-1- Accelerations 

For a one dimensional cylinder of length 2L and line density A, the net 
acceleration a distance d < L from the center of the cylinder, along the central axis 
of the cylinder, is 

: L+d GXdr f L ' d GXdr 

a>int{d) 

&int{d) 



/ jy«2 



' Jo 

L+d GX dr 



L-d 



r 



2 



a mt {d) = 7^-7^- ( 15 ) 



The net acceleration a distance d > L from the center of the cylinder, and 



13 



along the central axis of the cylinder, is 



a ex t{d) 
a e xt{d) 



L+d GX dr 

d-L r2 

GX GX 



d-L L + d 



(16) 



If the above described cylinder is perturbed such that its density becomes 
(1 + e) A within a distance of L\ of the center of the cylinder, then the acceleration 
a distance d from the center will be 



Qnet(^) 



GX _ _GA_ , GXe _ GXe A -C J 
L-d L+d 1 L x -d Li+d' « ^ -^1 



(17) 



GX GX _|_ GX t 



GXe 



L-d L+d 1 d-L 1 Li+d 



, L 1 < d < L. 



For the case of L\ < d < L, expanding to lowest order gives 



a nct (d) 
a ne t(d) 



GXi T 



'! + -) 



GX2df, eL 1 L 2 
1 + 



L 2 



rf 3 



(18) 



Figure [3] shows the accelerations of a perturbed cylinder as calculated from the 
above exact and approximate equations. For simplicity, A, G, and L\ are set to 1 
while e and L are set to 0.1 and 10 respectively. Just as with the disc, Equation 
[TBI accurately reproduces the accelerations at intermediate radii but underestimates 
the accelerations near L\ and L. 
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2.4-2. Timescales 



The square of the instantaneous collapse time is approximately 



The square of the global collapse time is thus 



C, = t\L), 



GX \ L 



K = 7n 1 + -T • (2°) 



For small perturbations 



— « h 

= lr (21 > 

Just as with the disc, the global collapse time is essentially the same in a perturbed 
and unperturbed cylinder and Equation [19] can be written as 

When e = 0, the local collapse time is essentially the global collapse time. The 
presence of a line density perturbation speeds up the local internal collapse. 

The square of the collapse timescale for a local region of size L\ is 



14 
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Equation [23] indicates that the smaller the perturbed region is in comparison to the 
total length of the cylinder, the faster the local collapse will be in comparison to 
global collapse. 

For the local collapse timescale to be significantly shorter than the global 
collapse timescale 

Li < JlL. (24) 

More specifically, for the local collapse timescale to be a factor of three smaller than 
the global collapse timescale, 

(I) 2 > * 

L x < 0.3V~eL. (25) 

If e = 0.1, then the perturbation length must be approximately 10 times smaller 
than the cylinder's length in order for the local collapse timescale to be three times 
smaller than the global collapse timescale. Taking L\ = Rj as the minimum size for 
a gravitationally unstable weak perturbation, this implies a minimum finite cylinder 
mass of a few tens of Jeans masses. A large perturbation, e = 1, would require 
a cylinder with approximately ten Jeans masses in order for the local collapse 
timescale to be three times smaller than the global collapse timescale. 

Figure H] shows the collapse timescales for a perturbed cylinder with 
Li — G — X — 1, e — 0.1, and L/L\ = 10. The solid line shows the timescales 



16 



calculated from the exact accelerations while the dashed line shows the timescales 
calculated from the approximate accelerations given by Equation [18j The collapse 
times for points interior to L\ are approximately a factor of three smaller than 
the collapse times for most points beyond L\. Just as for the disc, the collapse 
timescale abruptly drops to zero at L and Li due to the sharp jumps in density 
at these points. In more realistic filaments, these density changes would be more 
gradual and the actual collapse timescales would more closely resemble the collapse 
timescales derived from the approximate accelerations of Equation [18j As in the 2D 
case, the collapse timescale at L\ calculated from the approximate accelerations is 
roughly equal to the collapse timescales for most points interior to L±, as calculated 
from the exact accelerations. 



In the above models, the collapse timescales become infinitely small at density 
boundaries. This behavior is clearly unphysical and originates from our assumption 
that the disc and filament are infinitely thin. To investigate the significance of this 
assumption with respect to our results, we have calculated the accelerations along 
the central axis of a cylinder with uniform density p, total height 2L, and radius R. 
We find that the acceleration a distance d < L from the center of the cylinder, but 
still along the central axis, is 



2.5. Edge Effects - A Hybrid Collapse Mode 




(26) 



-17- 



Evaluating this at the edge of the cylinder, d = L, gives a finite acceleration of 



thus proving that it was the infinite thinness of our prior models that led to the 
infinitely small collapse timescales. 

As before, these accelerations can easily be turned into instantaneous collapse 
timescales. Figure |5] shows the collapse timescales for two cylinders, both with 
G — 1 and L = 10, where one cylinder has a radius R = 1 and uniform volume 
density p — 1 and the other cylinder is infinitely thin and has a line density 
equivalent to that of the other cylinder, namely A = ir. The collapse timescales of 
the two cylinders are essentially identical except for points very close to the end of 
the cylinders. The collapse timescales of the finite cylinder only exceed those of the 
infinitely thin cylinder by a factor of two or more for points within 0.02 L of the 
end. Since there is such a good agreement between the collapse timescales and since 
our global collapse timescales were calculated from our approximate accelerations, 
which did not trace the sharp increases in acceleration at the edges, our above 
analysis of the ratio of local to global collapse timescales should be valid for both 
infinitely thin and finite sized objects. 

Figure [5] also shows that while a finite radius cylinder does not have infinite 
accelerations at its edges, the timescale for collapse still decreases towards the edge 
in the outer part of the cylinder. For the particular set of parameters of the finite 
radius cylinder shown in Figure |5j the edge collapse timescale is roughly a factor of 
three less than the global collapse timescale, given by Equation |21j for a filament 
with this line density. This reduction in collapse time at the edge is due to the 




(27) 
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density contrast at the edge. 

A collapse timescale that decreases outwards leads to material at the edge 
running into material further inwards, which causes a build up of material. Such a 
density enhancement at the periphery of a cloud would then collapse on the local 
collapse timescale. We consider this edge collapse mode to be a third collapse 
mode that is a combination of a global and local collapse mode, since the initial 
accelerations are provided by the global gravitational potential while the final 
collapse only occurs within the local density enhancements along the edge of the 
cloud. 

The majority of a filamentary cloud should still collapse roughly on the global 
collapse timescale, rather than the edge collapse timescale, as the material inward 
of the edge will slow the edge's collapse. This is further supported by the fact that 
for the cylinders in Figure |5j the collapse timescales for over 85% of the cylinders' 
lengths are within a factor of two of the global collapse timescale. 

To test how important the sharp edges of our clouds are to the presence of 
an edge collapse mode, we examine a series of finite radius cylinders in which the 
density slowly decreases to zero, rather than abruptly dropping to zero, near the 
faces of each cylinder. We use a cylinder with R=l,L = 10,p=l, and G = 1 
as our comparison, sharp edged cylinder. To create a tapered cylinder, we set the 
central density of the cylinder, p , to be equal to 1 and outside of a distance z e d ge 
from the center of the cylinder, as measured along the central axis of the cylinder, 
we adopt a gaussian density profile: 




(28) 
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We examine five tapered cylinders where the values of a are 0.5, 1, 2, 4, and 6. 
Note that the case of a = 1 corresponds to the taper having a size similar to the 
radius of the cylinder. We extend the tapered edge out to a length of 3a, such that 
the density drops to less than 2% of the cylinder's central density. We choose the 
location where the tapered edge starts, z e d ge , by requiring that the total mass of 
the tapered cylinder is the same as that of our comparison sharp edged cylinder. 
We then numerically calculate the instantaneous accelerations along these cylinders 
and convert these accelerations into collapse timescales, as done previously in this 
paper. The density profiles and collapse timescales for these cylinders are shown in 
Figure [6j 

As the length scale of the taper increases, the collapse timescale towards the 
edge increases and the collapse timescale of the interior decreases, thereby reducing 
the effectiveness of the edge collapse mode, as expected. A decrease in collapse 
timescale towards the edge, however, is still very prominent in all of the tapered 
cylinders. The ratio of the minimum collapse timescale (edge collapse) to the 
interior collapse timescale (global collapse) is 0.32 for the sharp edged cylinder and 
only increases to 0.40 for the a = 1 cylinder. It is not until the a = 4 cylinder, in 
which the central, uniform density section of the cylinder is just 60% of the length 
of the sharp edged cylinder, before the ratio of the minimum collapse timescale to 
the interior collapse timescale rises above 1/2. We therefore suggest that this edge 
collapse mode should still be significant in realistically tapered filaments. While 
we have not performed a similar calculation for a finite thickness, tapered disc, we 
expect that a very large taper would also be required before this edge effect becomes 
negligible in a disc. 
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3. DISCUSSION 

3.1. Interpretation of Analytical Results 

We have shown that weak density perturbations in spheres, discs, and filaments 
will collapse faster than the entire cloud, but that the significance of these local 
collapse modes depends strongly on the geometry of the cloud. 

For spheres and discs, global collapse modes are significant. As shown in §2.21 
the timescale for the collapse of perturbations in spheres is only dependent upon 
the strength of the perturbation and thus, very large perturbations are required for 
a spherical cloud to fragment. The collapse timescale for perturbations in discs, 
as shown in §2.31 is dependent upon the relative length scale of the perturbation, 
in comparison to the total size of the disc, but since thermal support sets a 
minimum size scale for collapse, roughly a Jeans length, we find that discs must 
be unrealistically large, having masses greater than 10 4 Jeans masses, for weak 
perturbations to collapse reasonably quicker than the entire disc. As such, discs 
also require large perturbations, possibly seeded by supersonic turbulence, in order 
for significant fragmentation to occur. 

In §2.41 we find, however, that filamentary geometry provides a favorable 
situation for small perturbations to grow before global collapse overwhelms them. 
This situation occurs because the local collapse timescale in filaments has a stronger 
dependance on the length scale of the perturbation than in discs or spheres. A 
filament needs to be only a few Jeans lengths long in order for local collapse to 
occur much more rapidly than global collapse, even with only weak perturbations. 
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Thus, once filaments are formed, they are likely to fragment into multiple objects, 
independent of the nature of turbulence in the filaments, which suggests that star 
forming cores should readily occur in filamentary structures in molecular clouds. 
Recent observations have discovered a correlation between star forming cores and 
dense filaments, as discussed further in §3.2.21 

We have only examined the collapse timescales along the longest axis of our 
clouds and have assumed that our clouds are supported from collapsing along their 
short axis. This is not unreasonable as pressure support is most effective on small 
scales (i.e., scales smaller than the Jeans length). For an infinite, isothermal disc, 
regardless of the surface density, there is always an equilibrium configuration of 
material alon g the short ax i s in w hich the disc has a thickness of approximately a 



Jeans length (ISpitzer 



1942 



19781 ) and for isothermal spheres, there is no shorter 
axis and thus, no extra collapse mode to account for. For an isothermal, infinite 
filament, there is a critical mass per unit length for the filament to be stable against 



radial collapse and at this critical mass per unit 



of approximately a Jeans length ( IStodolkiewicz 



lengt 



1963 



i, the fi 



Ostriker 



ament has a radius 



1964). Below this 



Inutsuka fc Miyama 



19971 ). 



critical value, a filament can be pressure confined (e.g. 
but above this critical value, radial collapse cannot be prevented by thermal support 
alone. 

Our long axis collapse modes increase the mass per unit length of a filament 
and could easily raise the mass per unit length of a subcritical filament above the 
critical value. If such a filament contains a weak perturbation, the perturbation 
will reach the critical mass per unit length before the rest of the filament and thus, 
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radial collapse will occur first in the perturbation. This radial collapse will increase 
the speed at which the perturbation collapses and thus, will make fragmenting a 
filament even easier than what our analytic work suggests. 

While the radial collapse mode may determine the late time density evolution 
of a perturbation, collapse along the long axis should still determine the final mass 
of the object as the majority of the mass accessible to a perturbation is along the 
long axis of the filament. Since our analytic results show that perturbations on 
the smallest length scales collapse the quickest and since thermal pressure sets a 
minimum collapse length scale at the Jeans length, filaments should fragment to 
form roughly Jeans mass objects separated by a few Jeans lengths. 

This analytic work assumes that there are no effective large scale support 
mechanisms operating within a molecular cloud. If such a mechanism were 
present, molecular clouds could survive for multiple free fall times and even slow 
local collapse modes could occur. P ossible support mechanisms discussed in the 



literature include driy en turbulence (IBonazzola et al. 



Li & 


Nakamura 


2006 


) anc 


magi 


1976; 


Shu 


1977; 


Shu et al. 


1987) 



1992 



Krumholz et al 
n rr~ 



2006 



2006 J) and mag netic fields (IMouschovias fc Spitzerl Il976l ; iNakano 



Ballesteros-Paredes et al. 



( l2009af ) show that the galactic gravitational potential 
has a relatively negligible effect on a spherical molecular cloud, but for a filamentary 
cloud, the external poten t ial can be larger than the self-gravity of the cloud. 



Ballesteros-Paredes et al. 



( I2009bl ) examine the energy budget of the Taurus 
molecular cloud and find that over the entire cloud, the galactic tidal energy is 
larger than the gravitational energy and should act to disrupt the cloud. On the 
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smaller scales of clumps, however, 



Ballesteros-Paredes et al. 



fl2009bh find that the 



gravitational energy is much larger than the tidal energy and thus, in the Taurus 
molecular cloud, the galaxy's gravitational potential should prevent global collapse 
but should not significantly hinder small scale collapse. 



3.2. Comparison with Simulations and Observations 

In this section, we apply the results of our analysis to both simulations, in 
13.2. 1| and observations, in §3.2.21 



3.2.1. Simulations 
Simulations of supersonic turbulence show that the fragmentation of three- 



turbulence (e.g., 


Bate et al. 


2002. 


2003; 


Klessen et al. 


2000; 


Ostriker et al. 


2001 


Clark & Bonnell 


2004 




Li et al. 


2004 




Vazquez-Semadeni et al. 


2005 


). Magnetic 



fields are also only capable of providing support perpendicular to the field direction 
and thus, magnetically supported clouds will also tend to collapse to form sheets and 
filaments. As we have shown, once filamentary substructure forms, the filaments 
can readily fragment to form star forming cores. Thus, star formation may occur 
through a process of reduction of dimensions, wherein spherical and disc shaped 
clouds first fragment to form fil aments and t hen star formation occurs within these 
filamentary structures (see also lLarsonlll985l ). 



Sheet-like, or at least flattened, clouds would be expected if clouds assemble 
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in la rge-scale flows from the amb ient medium. Besides observational evidence (see 



e.g. 



Beaumont fc Williams 



20101 ) . a whole series of numerical experiments has 



demonstrated that the formation of clouds due to sweep-up by shells or flow collis ions 



can lead to flattened (yet 



Heitsch fc Hartmann 



2008 



turbulent) structures 



Banerjee et al. 



(jVazquez-Semadeni et al. 



2007 



20091 ) that, in so me cases, display ring- 



l ike st r uctures, as seen in the more idealized experiments by 



(200 



Burkert fc Hartmann 



Heitsch et al. 



( l2008af ) explicitly demonstrate that the ring-like filaments 
in their models are caused by global gravitational modes at the edge of the cloud. 
Yet from the current models it is difficult to see whether these ring structures will 
persist in a more general situation, where the colliding streams are not strongly 
spatially constrained. Unli ke the idealized, truly two-dimensional models of 



Burkert fc Hartmann! (120041 ) . the three-dimensional cloud- formation m odels do show 



local c ollapse modes induced by turbulence and/or strong cooling (see 



Heitsch et al. 



2008bl for a discussion of the timescales). 



The preferential formation of density enhancements at the edges of filaments 
and discs, as would be caused by the edge c ollapse mode described in £12.51 does 



occur in various numerical simulations (e.g., 



Bastien 



1983 



Hartmann fc Burkert 



20071 ) . In §2.51 we determined that while this edge collapse mode can be suppressed 



in filamentary structures by tapering the edges of the filament, the taper 



las to 



Li 



(2001) 



be quite large in order for the edge collapse mode not to be significant, 
examines the importance of sharp edges to the formation of a dense ring along the 
outer edge of a collapsing disc and finds that the formation of this ring is only 
suppressed when the edge of the disc is tapered and the length scale of the taper is 
at least comparable to the size of the uniform density portion of the disc. Similarly, 
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for filaments, the simulations of iNelson fc Papaloizoul (119931 ) show that prolate 
spheroids, which are highly tapered cylinders, do not always form fragments at their 
ends. 

Consideration of the short global free-fall time for finite molecular clouds 
suggests that it is hard to maintain a "quiescent" star formation mode whereby 
individual Jeans mass objects fall out of a "quasi-static" cloud. Rather, our results 
emphasize that the bulk of the cloud also partakes in the gravitational collapse. 
While observations of molecular clouds often refer to "quiescent" versus "clustered" 
star formation, the distinction in this work is between local and global collapse in 
a finite cloud. In our picture, global collapse inevitably occurs, but local collapse 
is more prone to happen at lower dimensions, namely in rings/edges for sheets, 
and at the ends of filaments. Yet, note that this kind of local collapse occurring 
at the edges of clouds is still at least partially driven by global accelerations in all 
cases. Likewise, local, pre-seeded perturbations are more likely to grow and collapse 
in clouds of lower dimensionality; perturbations in sheets, if large enough, and 
filaments have more of a chance to collapse before global collapse ensues. 

The similarity between the local and global collapse timescales in spheres and 
discs means that global collapse modes cannot be ignored in models or simulations 
of spherical or disc shaped molecular clouds. Thus, predictions of evolution 



timescales from models of cloud fragmentation in periodic b oxes ( Klessen et al. 



2000; 



2002 



Qstriker et al 



2001 



Heitsch et al. 



Vazquez-Semadeni et al 



2001 



Li et al 



2004 : 



Padoan fc Nordlund 



20051 ). or quasi-static approaches (IMyers 



2009), are 



strongly limited by their assumptions of boundary conditions and/or initial states. 



-26- 



Furthermore, since periodic box models do not incorporate physical processes 
occurring on scales larger than the box length, such as global collapse, it is 
exceedingly difficult to accurately run them for longer than the background global 
collapse timescale. As such, while periodic box simulations accurately follow the 
evolution of strong density perturbations, they cannot follow the collapse of small 
perturbations to completion. In other words, they only give meaningful results on 
timescales that are shorter than the global collapse timescale in the absence of the 
periodic boundaries. 



3.2.2. Observations 



Filamentary structure in star forming clouds appears to be ubiquitous. 
Filaments are present in both high mass and low mass star forming regions and 
are detected from the largest scales of molecular clouds all the way down to the 
scale of individual star forming cores; although, individual cores tend to be less 



elongated as thermal pressur e 



Schneider & Elmcgreen 



Kulkarni fc 



Jeiles 



Onishi et al. 



1996 



988 



1979 



plays a larger role on smalle r 



Scalo 



Loren 



Lada et al. 



1985 



1989 



1999 



Ballv et al. 



Nozawa et al. 



iiina et al. 



1987 



1991 



999 



spatial sea 



Hcv er et al. 



es e. 



1987: 



Tatematsu et a 



ieithausen et al. 



1993 



2002 



Hartmann 



2010 



2002 



Johnstone et al. 



Andre et al. 



2003 



Lada et al 



2007 



Mversl 12009 



Molinari et al 



20101 ). For instance, the Orion A molecular cloud is filamentary 



on large scales and contains filamentary substructure, inclu ding the integral 



shaped filament and the Orion Nebula Cluster (ONC) (e.g., 



Johnstone fc Bally 



Bally et al. 



1987 



19991 ). Such filamentary structure is seen in both gas tracers 
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and in star counts (e.g., iBally et al.lll987t iHillenbrand fc Hartmannlll998l ). Radial 
velocity measurements of stars and gas around the ONC have revealed a large 



scale velocity gradient, which has been interpreted as global co 



long axis of the filamentary structure of the ONC (IFuresz et al. 



lapse 



2008 



along the 



Tobin et al 



20091 ) . The large scale structure of the Orion A molecular cloud has also been 



remarkab 



y well reproduced by a sim ulation of a rotating disc with a density 



gradient (IHartmann fc Burkert 



20070 . 



The presence of filamentary structures appears to be closely tied to the 
star formation process, as young protostars and bound prestellar objects are 



preferentially locat ed within dense filaments (IHartmann 



Andre et al 



2002 



Andre et al 



2010|) 



(120 lOf ) further note that Class protostars and bound prestellar cores 
primarily exist in filaments that have mass per unit lengths larger than the critical 
mass per unit length at which an infinite, isothermal cylinder will be unstable to 
radial collapse (see §2. II) . suggesting that the radial collapse mode of filaments plays 
a vital role in fragmenting a molecular cloud. This association of star forming cores 
with filaments is what we would expect given that our analytical results indicate 
that only filaments can readily fragment without strong perturbations. 



4. CONCLUSIONS 

We have analytically estimated the global collapse timescales of spherical, disk 
shaped, and filamentary molecular clouds, as well as the local collapse timescales 
for small perturbations in these same clouds. We have shown that local collapse 
modes with timescales less than the timescales for global collapse do exist for 
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perturbed spheres, discs, and cylinders and that local collapse is more effective 
for lower dimensional objects. We find the square of the ratio of global to local 
collapse timescales is inversely dependent upon the square of the relative size of the 
perturbed region of a cylinder, is inversely dependent upon the relative size of the 
perturbed region of a disc and is independent of the relative size of the perturbed 
region of a sphere. Table [I] summarizes the most important derived equations 
relating to the collapse of perturbed spheres, discs, and cylinders. 

We find that filamentary geometry is the most favorable situation for the 
smallest perturbations to grow before global collapse overwhelms them and that 
filaments containing only a few Jeans masses and weak density perturbations can 
readily fragment. Conversely, we find that weak perturbations in realistically sized 
discs and spherical clouds do not collapse significantly faster than the entire cloud. 
Global collapse modes are significant in spheres and discs and cannot be ignored 
in simulations or models of such clouds. We also find that there exists an edge 
instability in discs and filaments that may lead to dense, locally collapsing structures 
along the periphery of these clouds. Our results suggest that star formation may 
occur through a process of reduction of dimensions, wherein spherical and disc 
shaped clouds first fragment to form filaments and then star formation occurs 
within these filamentary structures 
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Fig. 1. — Accelerations of a perturbed disc with R = 100, G — 1, E = 1, e = 0.1, 
and i?i = 1. The solid line shows the exact solution for the central accelerations 
while the dashed line shows the lowest order approximation to the exact solution. 
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Fig. 2. — Collapse timescales of a perturbed disc with i? = 100,G = l,S = l, 
e = 0.1, and R± — 1. The solid line shows the timescales calculated from the 
exact accelerations while the dashed line shows the timescales calculated from the 



approximate accelerations. 
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Fig. 3. — Accelerations of a perturbed cylinder with L\ — 1, G — 1, A = 1, e = 0.1, 
and L = 10. The solid line shows the exact solution for the central accelerations 
while the dashed line shows the lowest order approximation to the exact solution. 
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Fig. 4. — Collapse timescales of a perturbed cylinder with L 1 = 1,G = 1,A = 1, 
e = 0.1, and L = 10. The solid line shows the timescales calculated from the 
exact accelerations while the dashed line shows the timescales calculated from the 
approximate accelerations. 
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Fig. 5. — Collapse timescales of a cylinder with L — 10, G — 1, p — 1, and R = 1 
are shown as the solid line and the collapse timescales of an infinitely thin cylinder 
with L — 10, G — 1, and A = ir are shown as the dashed line. 
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Fig. 6. — Density profiles (left) and collapse timescales (right) of tapered cylinders. 
For all cylinders, G — 1, R — 1, and the central density is p = 1. The solid line is for 
a cylinder with L — 10 and sharp edges. The dotted, dashed, dot-dash, dash-triple 
dot, and long dashed lines correspond to cylinders where a is 0.5, 1, 2, 4, and 6, 
respectively, and where a is the standard deviation of the gaussian profile used for 
the tapered edge. The location of the tapered edge has been chosen so that all of the 
cylinders have the same total mass. Note how the collapse timescale still decreases 
significantly towards the edge, even for the most strongly tapered cylinder. 
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